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摘 要 : 准确 估算 草地 生物 量 ， 对 全 球 气 候 变 化 背景 下 的 陆地 生态 系统 碳 循环 研究 具有 重要 意义 .过 去 几 十 年 ， 
草地 生物 量 估算 研究 大 多 集中 在 北方 ， cd m 对 其 生物 量 进行 评估 的 
报道 较 少 。 本 文 以 云南 省 为 例 ， 应 用 2012—2014 年 草地 生物 量 野 外 调查 资料 和 同期 MODIS 3E REGIE, XE 
草地 地 上 生物 量 (AGB) 遥 感 估算 模型 ; 再 引入 草地 植被 群落 特征 (高 度 和 盖 度 ) 信 息 对 统计 模型 进行 优化 ， 并进 
行 生物 量 空间 反 演 。 结 果 表 明 : 优化 后 模型 的 估算 精度 由 原来 的 35.0% 提 升 为 43.7%; 反 演 得 到 云南 省 3 年 年 
均 AGB 的 总 量 介 于 1026.86 万 ~1 408.54 A t, 平均 为 1221.11 万 t; 从 空间 分 布 上 看 , 云南 省 草地 AGB 密度 
总 体 上 呈现 西部 高 东部 低 、 南 部 高 北部 低 的 格局 。 本 研究 首次 将 远 感 植被 指数 数据 与 实测 植被 群落 特征 参数 
结合 ,使 估算 精度 比 传统 的 纯粹 光学 遥感 模拟 方法 显著 提升 24.996, 但 精确 估算 大 面积 的 草地 AGB, 需要 进 
一 步 探 索 如 何 将 激光 雷达 数据 或 遇 感 立体 影像 中 提取 的 植被 特征 信息 应 用 于 草地 AGB 估算 研究 。 
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Estimation of grassland biomass using MODIS data and plant 
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Abstract: In the context of global climate change, the accurate estimation of grassland biomass is critical for terrestrial carbon 
cycling research. In China, most studies in this area have focused on grasslands in North China over the past decades. Only a 
few studies have estimated grassland biomass in South China, mainly due to difficulties in spatial complexity of plant species 
in the region. Therefore, it is necessary to develop a model for the estimation of grassland biomass in South China in order 
to analyze the spatial distribution of this vegetation type. In this study, we first developed a model for the estimation of 
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aboveground grassland biomass (AGB) in Yunnan Province using field sample and NDVI (normalized difference vegetation 
index) data (2012-2014), derived from MODIS sensor. The derived grassland characteristics (height and coverage) were then 
inputted into the model to improve the estimation accuracy. With the improved model, we used remote sensing and GIS 
platforms to map the spatial pattern of AGB in Yunnan Province. Finally, we carried out statistical analysis of AGB of 
grassland in a district in Yunnan Province and calculated the average density of AGB in multiple types of grassland. The 
results indicated that: 1) the model for the estimation of AGB of grassland was improved by the use of field data on plant 
community. Thus the goodness-of-fit (R°) of the model increased by 0.289 and the estimation accuracy of the model also 
increased (35.096—43.796) significantly. 2) During 2012-2014, annual total AGB in the study area was 1.03x107—1.41x10? tons, 
with an average value of 1.22x10? tons that accounted for 4.196 of total AGB in China. The results suggested that the area of 
grasslands in South China is not negligible. The density of AGB of grassland in Yunnan was highest in the eastern and southern 
regions of the province. 3) The density of AGB of grassland in the districts of Yunnan was 1 130.12-2 116.03 kg:hm ?. 
Grasslands with high AGB densities were in southern and southwestern areas of the province, including Xishuangbanna, 
Dehong and Puer. Grasslands with low densities were in northwestern and eastern areas of the province, including Diqing and 
Qujing. Moreover, AGB density of mutiple grassland types had a clear pattern, with an increasing trend from montane meadow 
to tropical herbosa. The order of the AGB density increase was: montane meadow (1 071.73 kg.hm ?) < lowland meadow 
(1552.45 kg-hm ?) < tropical shrub herbosa (1579.80 kg-hm ?) < warm-temperate shrub herbosa (1588.12 kg.hm?) < 
warm-temperate herbosa (1 771.02 kg-hm ?) < tropical herbosa (2 004.37 kg:hm ?). In the study, a remotely sensed vegetation 
index was first combined with field data on plant community. Using this approach, the accuracy of the results increased with 
24.9%, compared with the traditional approach which relies only on remote sensing data. Thus in order to improve the 
accuracy of AGB estimation for grasslands at a large scale, it was recommended that future studies attempt to incorporate 
grassland height derived from Light Detection and Ranging equipment (LIDAR) data or optical stereo images. 


Keywords: Yunnan Province; MODIS data; Aboveground biomass; Plant community characteristic; Space inversion 


草地 是 陆地 生态 系统 的 重要 组 成 部 分 ， 在 全 球 
气候 变化 和 碳 循环 中 扮演 着 重要 角色 趾 。 草 地 生物 
量 可 以 用 来 判断 草地 的 生长 状况 和 生产 潜力 ， 准 确 


TARRE, @ 基 于 生物 量 扩 张 /转换 因子 或 系数 ; 
外 定点 数据 与 下 感 数据 结合 。 在 以 上 方法 中 , E 
点 数据 与 遥感 数据 结合 的 方法 在 大 面积 估算 植被 生 


估算 草地 生物 量 及 其 动态 变化 对 研究 全 球 碳 循环 和 
评估 草地 植被 的 生态 效益 具有 重要 意义 1。 

我 国 草地 面积 为 3.55x10* hm2， 占 国土 面积 的 
41.794", Bj 20 世纪 90 年 代 中 期 , 一 些 学 者 相继 开 
展 了 我 国 草地 植被 生物 量 估算 的 研究 6。 这 些 研 究 
主要 集中 在 北方 区 域 中 ,包括 西藏 、 内 蒙古 、 青 海 及 
新 疆 等 地 。 如 牛 志春 等 "利用 TM 遥感 图 像 和 同期 
野外 实测 数据 ,建立 了 青海 湖 环 湖 地 区 草地 植被 生 
物 量 的 遥感 监测 模型 ; 张 峰 等 nM 利用 内 蒙古 典型 草 
原 连 续 13 年 的 地 上 生物 量 (aboveground biomass, 
AGBA, 基于 遥感 信息 进行 生态 系统 碳 循环 过 程 
的 CASA 模型 验证 ， 并 分 析 了 内 蒙古 典型 草原 1982— 
2002 年 植被 净 初 级 生产 力 的 时 间 变 异 及 其 影响 因 
子 。 据 统计 ， 南 方 各 省 有 草山 草 坡 6000 多 万 hm, 
占 全 国 草地 面积 近 1/600， 在 北方 草地 趋 于 退化 的 
情形 下 PM” 3 南方 草地 在 全 国 草地 生态 系统 中 的 地 
位 越发 重要 。 然 而 ， 南 方 地 形 复杂 、 气 候 差 异 显著 
等 因素 导致 草地 的 空间 异 质 性 非常 明显 ， 表 现 为 草 
地 类 型 繁多 ， 且 分 布 零散 ， 因 而 对 云南 省 草地 生物 
量 进行 反 演 与 估算 的 研究 尚 属 空白 5。 

生物 量 估算 模型 一 般 为 : 四 基于 过 程 模型 ; OH 


物 量 时 应 用 较 多 , 尤其 在 具有 野外 生物 量 调查 数据 
库 的 区 域 。 现 有 的 清单 数据 和 野外 调查 数据 在 进 
行 植被 生态 系统 年 度 动态 评估 时 存在 缺陷 , 包括 缺 
eaa 间 覆 盖 、 较 低 的 空间 分 辩 率 、 行 政 边 界 
间 不 连续 和 明显 的 时 间 党 后 性 等 " 。 明 感 数据 具有 
广泛 的 空间 覆盖 、 较 高 的 时 空 分 辩 率 及 良好 的 一 致 
性 和 稳定 性 " ”可 以 从 一 定 程度 上 弥补 上 述 不 
足 。 植 被 遥感 定量 研究 利用 已 有 的 实地 调查 数据 ， 
建立 以 遥感 植被 指数 或 者 环境 因子 为 自 变 量 、 生 物 
量 为 因 变 量 的 回归 模型 ， 在 此 基础 上 推算 整个 区 域 
BEHE., 在 利用 遥感 影像 (如 Landsat 和 MODIS 
数据 等 ) 计 算 的 各 种 植被 指数 中 ， 归 一 化 植被 指数 
(normalized difference vegetation index, NDVD) 是 较 
为 常用 的 植被 生物 量 指示 指标 , 为 监测 不 同 尺 度 植 
被 生物 量 的 时 空 分 异 特 征 提供 了 一 个 理想 的 可 量化 
参数 请”]。 但 是 , 用 NDVI 估算 生物 量 有 一 定 的 局 
限 性 , 包括 遥感 平台 及 传感器 性 能 限制 、 影 像 容易 
受到 云 和 气 溶胶 的 影响 以 及 NDVI 存在 数据 饱和 的 
问题 等 %“。 前 人 研究 发 现 ， 生物 量 随 着 物种 组 成 和 
植被 高 度 的 变化 而 变化 1， 激光 雷 达 具 有 直接 测量 冠 
层 高 度 分 布 的 特点 ,适合 大 面积 生物 量 估算 PM" 因此 ， 


http://www.ecoagri.ac.cn 


532 


中 国生 态 农业 学 报 2017 58-2578 


将 植被 群落 特性 信息 (如 ， 高 度 和 盖 度 等 ) 与 NDVI 联 
合 ， 可 以 为 植被 生物 量 研 究 提 供 新 的 思路 。 本 研究 
考虑 到 雷达 数据 的 可 获取 性 以 及 植被 高 度 跟 种 群 的 
相关 性 5 在 进行 草地 生物 量 反 演 时 加 入 草地 植被 
群落 特征 参数 ,为 用 模型 进行 大 面积 生物 量 估算 提 


供 参 考 。 


本 文 以 云南 省 为 研究 区 域 , 利用 MODIS 数据 计 
同期 草地 AGB 实测 数据 ， 建 立 基 


PE PAN 


算 的 NDVI, 结合 
T NDVI 的 草地 AGB 估算 模型 ， 并 应 用 可 


高 度 和 盖 度 等 特征 ， 对 估算 模型 进行 优化 ,探讨 d 
将 植被 群落 特征 信息 引入 草地 AGB 估算 模型 是 否 
有 助 于 提升 估算 精度 ; ORMA Fit AGB 的 空间 分 
布 规律 ; 以 云南 为 代表 的 南方 草地 在 我 国 草地 生 
态 系统 中 的 地 位 。 旨 在 为 我 国 草地 生物 量 估算 、 载 


冀 能 力 评估 以 及 持续 利用 提供 科学 依据 。 


1 数据 获取 与 研究 方法 
11 ”研究 区 概况 

云南 省 地 处 中 国 西南 边 隆 , 位 于 21?09'~29?15'N、 
97?"30'~106?*00'E( 图 1)， 区 域 海 拔 差 异 显 著 ， 为 76.4~ 
6 740 mo 云南 省 总 土地 面积 为 3.94x10" hm’, 其 中 
草地 面积 为 1.527x10" hm”， 占 总 土地 面积 的 38.7%, 
可 利用 草地 面积 为 1.187x10 hm, SRX E F 
高 原型 的 热带 季风 气息， 具有 年 温差 小 , 日 温差 大 ， 
雨量 充沛 的 气候 特征 。 全 省 平均 气温 24 C, FP 
降雨 量 1 100 mm, 冬 春季 降水 仅 占 年 降水 的 20% 左 
t. 良好 的 自然 环境 为 植物 生长 提供 了 适宜 的 光照 、 
土壤 、 气 候 等 环境 ,形成 了 生物 资源 的 多 样 性 。 根 
据 实地 调研 情况 ,云南 省 的 植被 类 型 图 主要 有 高 塞 
草 旬 类 、 温 性 草原 类 、 山 地 草 旬 类 、 暖 性 ( 灌 ) 草 从 类 、 
热 性 ( 灌 ) 草 从 类 、 干 热 稀 树 灌 草 从 类 等 几 大 类 。 


直 物 群 落 的 
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1 研究 区 草地 生物 量 地 面 样 点 分 布 图 


Fig. 1 


12 ”研究 方法 
1.2.1 草地 生物 量 地 面 调 查 

外 业 调 查 时 间 主 要 集中 在 2012 一 201 
KZ, 根据 草地 类 型 的 空间 分 布 特征 和 


Distribution of sampling site 


s of investigation of aboveground grassland biomass in the study area 


共 在 云南 地 区 设计 了 305 块 调 查 样 地 。 在 2012— 
2014 年 8 一 10 月 期 间 ， 对 分 布 面积 较 大 的 草地 进行 
了 调查 , 采样 点 包括 云南 省 58 个 县 /县 级 市 。 在 进行 
数据 处 理 时 剔除 了 表现 异常 的 个 别 采 样 点 数据 ， 最 


4 年 草地 生 
面积 大 小 ， 
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终 剩 下 684 个 数据 用 于 建 模 和 验证 (图 1)， 其 中 2012 
年 设置 138 个 采样 点 ; 2013 年 设置 230 个 样 点 ; 2014 
年 扩大 了 草地 调查 的 范围 和 类 型 ,设置 了 316 TX 
样 点 。 

样 地 设置 在 草地 植被 空间 分 布 均一 ， 具有 代表 
性 的 典型 地 段 。 在 植被 群落 均一 性 较 好 的 样 地 内 设 
iB 1 或 2 个 0.5 mx0.5 m 的 样 方 , 在 植被 类 型 比较 复 
杂 、 分 布 不 均匀 的 区 域 设 置 3 或 4 个 1.0 mx1.0m 的 
样 方 。 在 样 地 的 基本 特征 调查 中 ,主要 记录 样 地 所 
隶属 行政 区 、 草 地 类 型 、 地 形 、 季 节 利 用 方式 和 利 
用 状况 等 。 利 用 GPS 测定 样 方 的 经 纬度 和 海拔 ， 同 
时 在 样 方 内 采用 常规 植被 调查 法 测定 植物 种 数 、 植 
被 盖 度 、 群 落 平 均 高 度 及 总 地 上 生物 量 等 指标 。 
1.2.2 ”遥感 数据 获取 及 处 理 

本 研究 所 使 用 的 遥感 数据 为 USGS( 美 国 地 质 调 
查 局 ) 提 供 的 EOS/MODIS MOD13Q1 数据 ， 该 数据 
空间 分 辨 率 为 230 m, 是 16 d 合 成 图 像 , 数据 格式 为 
EOS-HDF， 空 间 位 置 在 全 球 正弦 曲线 投影 SIN(sinusoidal 
projection) 系 统 中 编号 为 h26v06 和 h27v06 的 两 景 图 
R, 时间 序列 为 2012—2014 年 草原 生长 季 (8 一 10 
H), 共计 6 幅 云 量 小 于 10% 的 遥感 影像 。 

利用 MODIS 专业 处 理 软件 MRT(MODIS Reprojection 
Tools) 对 下 载 的 NDVI 数据 进行 投影 和 格式 转换 处 
理 ， 并 设置 为 统一 的 Albers 投影 (第 1 标准 纬度 : 
25?N, 第 2 标准 纬度 : 47*N， 中 央 经 线 : 105?E, KH 
基准 面 : WGS84)。 根 据 研究 区 植被 特征 ， 并 参考 国 
内 外 草地 地 上 生物 量 相 关 研 究 ， 利 用 遥感 软件 
ENVI 及 地 理 信息 系统 软件 ArcGIS, M MODIS 数据 
中 提取 归 一 化 植被 指数 NDVI, 计算 公式 ”为 : 


NDVI = PNR ^ PR (1) 
PNIR + PR 


xU: pum 为 近 红外 波段 的 反射 率 ，A 为 红 光 波段 
的 反射 率 。 

植被 覆盖 度 与 NDVI 之 间 具 有 极 显著 线性 相关 关 
系 , 因此 可 以 利用 NDVI 提取 植被 覆盖 度 信息 Fl。 
根据 像 元 二 分 模型 ,每 个 像 元 的 NDVI 值 可 以 表达 
为 由 绿色 植物 部 分 所 贡献 的 信息 NDVI,, 与 裸 土 部 
分 所 贡献 的 信息 NDVI,2 部 分 组 成 ， 则 植被 覆盖 度 
的 公式 中 如下: 

NDVI- NDVI,.; 
f= NDVI e - NDVI {oi Qe 

式 中 : fo 为 植被 覆盖 度 ; NDVI 为 完全 裸 土 或 无 杆 
被 覆盖 区 域 的 NDVI 值 ， 本 研究 中 选取 土地 单元 内 累 
只 频率 为 0.5% 的 NDVI 值 为 NDVIsor; NDVIvee 则 表示 
完全 被 植被 覆盖 区 域 的 NDVI 值 ， 本 研究 中 选取 土地 


单元 内 累积 频率 为 99.5% 的 NDVI (F7 NDVIvego 
1.2.3 ”回归 模型 建立 与 精度 验证 

利用 统计 软件 SPSS 进行 NDVI、 植 被 特征 信息 
(高 度 和 盖 度 ) 与 样 点 实测 AGB 数据 的 相关 分 析 和 回 
归 分 析 , 得 到 以 NDVI 植被 特征 为 自 变量 的 AGB XÉ 
感 估算 回归 模型 。 然 后 , 利用 实测 样 点 AGB 数据 与 
建立 的 遥感 估算 模型 计算 的 AGB 进行 比较 , 通过 均 
方 根 误差 RMSE) 和 估算 精度 (Accuracy) 对 模型 进行 
HENDI! RMSE 和 Accuracy 的 计算 公式 如 下 : 


XT- 
RMSE = (3) 


RMSE 
Y 
式 中 : RMSE 为 均 方 根 误差 ; Accuracy 7Jfh T8, Y, 
为 实测 草地 AGB(gm-), 7 为 草地 AGB 估算 值 (gm- )， 

Y 为 实测 草地 AGB 平均 值 (g-m”), N 为 样 点 数 。 


2 结果 与 分 析 
2.1 AGB-NDVI 建 模 与 检验 

从 2012 一 2014 年 地 面 样 方 调查 的 684 个 AGB 数 
据 里 面 随机 挑选 出 600 个 数据 , 用 ArcGIS 从 MODIS 
数据 中 提取 与 地 面 样 方位 置 对 应 的 同时 期 NDVI 值 ， 
建立 草地 单位 面积 的 AGB 与 NDVI 的 一 元 线性 回归 
方程 。 结 果 发 现 AGB 与 NDVI 呈 显 著 正 线性 相关 
(R=0.297, P«0.01), 回归 方程 为 : 

y=3 331.172x-405.647 (5) 

式 中 :y 为 AGB,x 为 NDVI。 

将 预 留 的 84 个 地 面 样 方 数据 代入 公式 (3) 和 (4) 
对 回归 方程 (5) 进 行 精度 验证 ,发 现 估 算 精 度 为 35.0% 
(图 2)。 


Accuracy = (1— )x10096 (4) 
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实测 值 Measured value (kg:hm?) 


2 云南 省 草地 地 上 生物 量 (AGB) 基 于 NDVI 的 估算 精 
度 验 证 
Fig.2 Accuracy validation for grassland aboveground 
biomass estimation based on NDVI in Yunnan Province 
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22 不 同类 型 草地 植被 特征 与 生物 量 

据 从 云南 省 获得 的 草地 AGB 野外 调查 资料 可 
知 ， 云 南 省 主要 的 草地 植被 类 型 有 9 大 类 (图 3)。 从 
3 可 以 看 出 ， 随 着 草地 类 型 由 高 寒 到 温 性 、 再 到 
热 性 的 变化 过 程 中 ， 草 地 植被 的 平均 高 度 (average 
height，AH) 逐 渐 增 加 ,平均 盖 度 (average coverage, 
AC) 没 有 明显 变化 ,而 草地 植被 的 AGB. 呈现 出 逐渐 
增 大 的 趋势 。 此 外 , 草地 AGB 和 AH 的 标准 误差 较 
大 , 而 AC 的 标准 误差 较 小 ,说 明 在 同类 草地 植被 中 ， 
AGB 和 AH 的 变化 波动 较 大 , AC 的 变化 较 小 。 进 一 
步 对 AGB 与 AH 和 AC 做 相关 性 分 析 , 结果 显示 
AGB & AH, AC 均 有 明显 的 正 相关 关系 ( 表 1)。 


度 Height (cm) 


高 


AGB (kg:hm?) 
E 
© 


1 2 3 4 5 6 7 8 9 
草地 类 型 Grassland type 


3 云南 省 不 同 草地 植被 类 型 的 高 度 、 盖 度 及 地 上 部 生 
物 量 (AGB) 
Fig.3 Height, coverage and aboveground biomass (AGB) of 
different types of grassland in Yunnan Province 
1: 高 寒 草 多 类 ; 2: 高 寒 草 旬 草原 类 ; 3: 山地 草 旬 类 ; 4: HE 
性 草丛 类 ; 5: 暖 性 灌 草丛 类 ; 6: 温 性 草 旬 草原 类 ; 7: 热 性 灌 草 
从 类 ; 8: 热 性 草丛 类 ; 9: 干 热 稀 树 灌 草 从 类 。1: alpine meadow 


type; 2: high-cold meadow-steppe type; 3: montane meadow type; 4: 


warm-temperate herbosa type; 5: warm-temperate shrub herbosa 
type; 6: temperate meadow-steppe type; 7: tropical shrub herbosa 
type; 8: tropical herbosa type; 9: arid-tropical shrub herbosa 
scattered with tree type. 


R1 云南 省 草地 地 上 部 生物 量 (AGB) 与 植被 高 度 (AH) 
和 盖 度 (AC) 的 相关 分 析 
Table 1 Correlation analysis between aboveground biomass 
(AGB) and height (AH) and coverage (AC) of grassland in 
Yunnan Province 


AGB-AH AGB-AC 
相关 系数 . 0.539 0.104 
Correlation coefficient (R) 
P 0.000 0.004 


此 ， 融 入 草地 植被 的 群落 特征 (高 度 和 盖 度 ) 有 助 于 
更 准确 地 估算 草地 AGB。 
2.3 草地 AGB 估算 模型 优化 

将 植被 群落 特征 信息 加 入 回归 模型 , 用 600 个 
样 点 进行 以 AGB 为 因 变量 , 以 NDVI、 高 度 和 盖 度 
为 自 变 量 的 多 元 线性 回归 分 析 ， 结 果 发 现 回 归 方 程 
的 拟 合 优 度 有 了 显著 提升 (R=0.614，P<0.01),， 通过 
检验 发 现 各 自 变 量 之 间 不 存在 多 重 共 线 性 ( 表 2)。 回 
归 方 程 为 : 

y=2 280.171x1+36.842x2+15.892x3-2 015.156 (6) 

式 中 : y X AGB, xi X NDVI, x; 为 植被 平均 高 度 ，x3 
为 植被 盖 度 。 

将 预 留 的 84 个 地 面 样 点 数据 代入 公式 (3) 和 (4) 
对 回归 方程 (9) 进 行 精度 验证 ， 发 现 估算 精度 提升 为 
43.7%( 图 4)。 优 化 后 的 回归 模型 比 之 前 的 拟 合 优 度 
提升 了 0.289， 估 算 精 度 提高 了 8.7%， 显 然 ， 模 型 的 
优化 对 于 准确 估算 云南 省 草地 AGB 具有 很 好 的 促 
进 作 用 。 


R2 云南 省 草地 地 上 部 生物 量 估算 回归 模型 系数 及 其 
检验 结 
Table 2 Coefficient of equation and its test result for 
estimation of grassland aboveground biomass in Yunnan 


Province 
系数 D 方差 膨胀 
Coefficient 因子 VIF 
归 一 TE] A 
归 一 化 植被 指数 2 280.171 6.107 0.000 1.057 
NDVI 
mE Height 36.842 16.178 0.000 1.051 
盖 度 Coverage 15.892 4.883 0.000 1.019 
常数 项 
常数 项 -2 015.156 -5.078 0.000 = 
Constant term 
24 BAE AGB 反 演 与 空间 分 布 特征 


在 利用 优化 后 的 回归 模型 对 云南 省 草地 AGB 
进行 反 演 时 ,考虑 到 全 局 空间 范围 草地 植被 群落 特 
性 数据 的 可 获取 性 ， 在 本 研究 中 进行 了 数据 的 简化 
处 理 ， 即 用 不 同类 型 草地 植被 的 平均 高 度 代替 真实 
高 度 ， 然 后 代入 回归 模型 应 用 GIS 中 的 栅 格 计算 器 
等 技术 手段 进行 空间 反 演 。 
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4 云南 省 草地 地 上 生物 量 (AGB) 基 于 NDVI 和 植物 特 
性 的 估算 精度 验证 
Fig.4 Accuracy validation for grass aboveground biomass 
estimation based on NDVI and plant characteristics (height and 
coverage) in Yunnan Province 


1£ 1 : 100 万 草地 类 型 图 分 类 系统 中 ， 未 对 云南 
省 的 高 寒 草 甸 、 高 寒 草 旬 草 原 、 温 性 草 旬 草原 和 干 
热 稀 树 灌 草丛 这 4 类 较 少 的 草地 类 型 进行 统计 ， 只 


表 3 云南 省 主要 草 


分 为 6 种 类 型 : 低地 草 旬 类 、 山 地 草 旬 类 、 暖 性 草丛 
类 、 暖 性 灌 草丛 类 、 热 性 草丛 类 及 热 性 灌 草 从 类。 其 
中 , 除了 低地 草 旬 类 外 ,其 余 类 型 草地 的 高 度 在 地 面 
样 方 调查 中 均 有 统计 (图 3); 根据 以 前 的 研究 5 0， 低 
地 草 旬 的 高 度 为 18.77~31.00 cm， 取 其 平均 值 ， 得 到 
低地 草 旬 的 AH 为 24.89 cmo 最终 得 到 云南 省 6 种 
主要 草地 类 型 的 AH 及 其 标准 误 ( 表 3)。 

根据 从 MODIS 数据 中 提取 的 2012 一 2014 年 的 
NDVI 数 据 , 利用 ArcGIS 计算 3 年 NDVI 的 最 大 值 、 
最 小 值 和 平均 值 ( 图 5); AH 的 最 大 值 为 均值 减 去 标 
ÆR, AH 的 最 小 值 为 平均 值 减 去 标准 误 ; 根据 
NDVI, 可 以 计算 得 出 植被 覆盖 度 的 最 大 值 、 平 均值 
和 最 小 值 (图 0)。 再 应 用 GIS 中 的 栅 格 计算 器 将 
NDVI 值 、AH 及 植被 覆盖 度 代 入 回归 模型 (60)， 进 行 
云南 省 草地 生物 量 的 空间 反 演 ， 得 到 3 年 AGB 密度 
的 最 大 值 、 最 小 值 和 平均 值 (图 7)。 可 以 看 出 , 云南 
省 草地 AGB 具有 明显 的 空间 分 布 特征 , 西部 AGB 
密度 普 逗 比 东 部 高 ,南部 总 体 上 比 北 部 高 ， 这 可 能 
与 气温 、 降 水 等 环境 因素 的 影响 有 关 。 


地 类 型 的 植被 特征 


Table 3 Vegetation characteristics of main grassland types in Yunnan Province 


平均 高 度 Average height (cm) 平均 盖 度 Average cover (96) 


草地 植被 类 型 
Grassland type 均值 标准 误差 均值 标准 误差 
Mean value Standard error Mean value Standard error 

低地 草 旬 类 Lowland meadow type 24.89 一 78.73 一 
山地 草 旬 类 Montane meadow type 12.60 0.70 92.01 0.80 
暖 性 草 从 类 Warm-temperate herbosa type 21.30 2.04 71.67 3.43 
暖 性 灌 草丛 类 Warm-temperate shrub herbosa type 23.78 2.99 82.44 4.71 
热 性 草丛 类 Tropical herbosa type 33.21 1.26 84.07 0.92 
热 性 灌 草丛 类 Tropical shrub herbosa type 29.62 1.78 84.94 0.91 


云南 省 2012—2014 年 草地 AGB 密度 介 于 1 376.76~ 
1 888.51 kghm”, 平均 为 1 637.21 kg:hm“。 云 南 省 
总 的 草地 面积 为 7 458 498.06 hm”， 由 此 计算 得 到 ， 
年 均 总 的 草地 AGB AF 1 026.86 万 ~1 408.54 F t, 
平均 为 1 221.11 万 t。 从 行政 分 区 上 看 ( 表 4), 密度 
较 大 的 草地 位 于 西双版纳 、 德 宏和 普洱 等 地 ， 主 要 
分 布 在 云南 的 南部 和 西南 部 地 区 ; 而 密度 较 小 的 草 
地 位 于 西北 的 迪 庆 和 东部 的 曲靖 等 地 。 总 AGB 最 大 
的 地 区 位 于 红河 哈尼 族 莱 族 自治 州 , 达 148.16 H t, 
其 次 是 普洱 、 昭 通 等 地 ， 而 总 AGB 最 小 的 地 区 位 于 
德 宏和 怒江 傈 傈 族 自治 州 等 地 。 

经 过 统计 发 现 , 不 同 草地 类 型 的 AGB 差异 较 大 
(R 5)， 从 图 1 的 草地 类 型 分 布 图 中 可 以 看 出 ，AGB 
密度 最 小 的 山地 草 旬 类 主要 分 布 在 云南 的 西北 部 ， 


而 最 大 的 热 性 草丛 和 暖 性 草丛 主要 分 布 在 中 部 和 南 
部 。 暖 性 草丛 和 暖 性 灌 草丛 分 布 最 为 广泛 ,这 两 类 
草地 的 总 体 AGB BRA, 分 别 达 611.26 万 t 和 
412.92 万 to 
3 Wit 

TARR se S36 ^E 5 38 BS — DERRE, 
不 同类 型 的 草地 植被 之 间 ， 高 度 和 盖 度 等 植被 特征 
存在 明显 的 差异 。 徐 敏 云 等 9 对 河北 省 草地 盖 度 、 
高 度 和 产 草 量 的 年 际 变化 进行 了 监测 研究 ,发现 7 
种 类 型 的 草地 植被 之 间 的 盖 度 差异 很 大 ， 其 中 山地 
草 旬 的 植被 盖 度 最 大 ,为 93.97%; 温 性 草原 的 植被 
盖 度 最 小 , 为 57.63%。 草 地 之 间 的 高 度 差 异 也 比较 
显著 ， 暖 性 草丛 草地 植被 高 度 最 高 ， 为 4925 cm; 
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NDVI of grassland from 2012 to 2014 in Yunnan Province Province 


低地 草 旬 最 低 , 为 18.77 cm。 HEU 根据 气候 - 植 。 形 地 貌 最 复杂 多 变 的 地 区 之 一 ; 全 省 海拔 相差 很 大 ， 
被 特征 ,选取 森林 草原 、 典 型 草原 和 荒漠 草原 3 种 ë 最 高 点 与 最 低 点 相差 6 663.6 mW。 特殊 的 地 形 地 貌 
不 同类 型 草地 进行 群落 特征 研究 ,结果 也 表明 3 种 和 悬殊 的 海拔 差异 形成 了 云南 多 样 性 气候 ,控制 了 
草地 的 盖 度 和 高 度 存 在 显著 性 差异 。 云 南 省 山地 面 。 复杂 的 水 热 条 件 引 起 的 物种 选择 、 资 源 竞 争 和 生境 
只 约 占 94%， 且 地 形 陡 肉 、 切 割 破碎 中 是 全 球 地 WAL, 造成 了 生态 系统 的 物种 多 样 性 5?， 由 此 形成 
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大 值 (a)、 平 均值 (b) 和 最 小 值 (¢) 
Fig.7 Maximum (a), mean (b) and minimum (c) values of 
aboveground biomass (AGB) of grassland from 2012 to 2014 in 
Yunnan Province 


草地 群落 间 显 著 的 空间 异 质 性 。 本 研究 ( 表 3) BB, 
云南 不 同类 型 的 草地 植被 盖 度 和 高 度 差异 明显 , 与 
徐 敏 云 等 n4 的 结果 基本 一 致 。 

草地 植被 的 高 度 、 盖 度 等 特征 是 描述 群落 生长 
状况 、 反 映 产量 高 低 较为 理想 的 特征 量 *I。 大 多 数 


表 4 云南 省 各 市 / 州 2012 一 2014 年 草地 平均 地 上 部 生物 
量 (ABG) 
Table 4 Averages of above ground biomass (AGB) of 
grassland of cites in Yunnan Province from 2012 to 2014 


市 (自治 州 ) 草地 面积 AGB 总 AGB 
E coa e NEC rn 
昆明 Kunming 356 902.24 439.78 51.39 
曲靖 Qujing 584 349.21 343.79 78.52 

EZ Yuxi 327 826.73 663.43 54.53 
昭通 Zhaotong 588 218.43 530.57 90.03 
楚雄 Weichu 471 578.27 531.66 72.23 
红河 Honghe 877 068.51 689.25 148.16 
文山 Wenshan 421 052.99 605.48 67.60 
普洱 Puer 722 718.19 998.01 144.40 
BE 410 845.97 2 116.03 86.94 
X3 Dali 491 058.55 593.07 78.23 
保山 Baoshan 318 075.86 826.86 58.11 
1&7: Dehong 182 384.74 2 033.23 37.08 
丽江 Lijiang 403 010.93 519.37 61.23 
怒江 Nujiang 263 080.30 623.76 42.72 
迪 庆 Diqing 608 272.92 130.12 68.74 
临 沧 Lincang 432 054.24 881.79 81.30 
总 计 Total 7 458 498.06 — 1 221.21 


X5 云南 省 各 市 / 州 2012 一 2014 年 不 同类 型 草地 平均 地 
上 部 生物 量 (AGB) 
Table 5 Averages of above ground biomass (AGB) of different 
types of grassland in Yunnan Province from 2012 to 2014 


H 只 总 
eea NE NEN 
rassland area -2 ota 
kg:h 
Grassland type (hm?) (kg-hm ) (100) 
j E" 
低地 草 甸 类 667354.90 | 1552.45 103.60 
Lowland meadow 
g 类 
山地 草 旬 类 521 531.62 1071.73 55.89 
Montane meadow 
i 从 类 
暖 性 草丛 类 3 451 442.61 1 771.02 611.26 
Warm-temperate herbosa 
暖 性 灌 草丛 类 
Warm-temperate shrub 2600033.82 1588.12 412.92 
herbosa 
H 从 类 
热 性 草丛 类 72 091.93 2004.37 14.45 
Tropical herbosa 
HUME EE M 
热 性 灌 革 从 类 276 218.11 1 579.80 43.64 
Tropical shrub herbosa 
总 计 Total 7 458 498.06 一 1 221.21 


草地 植被 高 度 增 长 与 生物 量 增 长 呈正 相关 ， 生 物 量 
峰值 与 植株 高 度 峰 值 同 时 出 现 中 。 本 研究 中 植被 高 
度 与 生物 量 呈 显著 的 正 相关 (R=0.539,，P<0.001)。 草 
地 植被 盖 度 也 是 生物 量 高 低 的 重要 表征 量 ， 陈 生 云 
SE FERREA NÉ ARCU Fb b RETE RU AE SE B X REN, A 
现 群 落 盖 度 与 生物 量具 有 极 显著 的 线性 相关 关系 
(R-0.96, P«0.001), 与 本 研究 结果 相符 (R=0.104，P< 
0.001, X 2)。 因此 ， 本 研究 综合 MODIS 数据 与 实测 
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的 不 同类 型 草地 特征 数据 进行 草地 AGB 的 估算 , 使 
估算 精度 比 传统 的 光学 遥感 模拟 方法 显著 提升 了 
24.9%， 而 此 方法 对 于 准确 估算 草地 AGB 的 应 用 尚 
少见 报道 。 大 量 研 究 证 实 森 林 垂 直 结 构 参 数 的 反 演 
有 助 于 森林 生物 量 的 估 测 ， 于 上 颖 等 个 通过 对 大 光斑 
激光 雷达 数据 GLAS 波形 的 处 理 ， 建 立 模型 反 演 森 
林 生 物 量 ,理论 精度 达到 91.3%。 本 研究 方法 为 今后 
估算 草地 生物 量 提供 了 一 个 新 的 思路 。 

结果 显示 云南 省 2012—2014 年 草地 AGB 均值 
为 1221.11 万 t, 在 相同 的 草地 面积 下 ,， 朴 世 龙 等 名 
在 全 国 尺 度 估 算 的 云南 省 草地 生物 量 比 本 研究 估算 
的 均值 低 185.11 万 t; 在 区 域 水 平 上 ，Xu 等 (和 1 利用 
MODIS 影像 和 地 面 实测 数据 估算 了 全 中 国 6 大 草地 
区 域 的 生产 力 ， 得 出 云南 省 2004 的 估算 结果 为 
1 520.33 万 1,2005 年 为 1 154.87 5 t, 与 本 研究 估算 
结果 接近 。 由 于 Xu 等 采用 的 云南 省 草地 面积 
(104 274 km 大 于 本 研究 的 草地 面积 ， 因 此 计算 的 
单位 面积 生物 量 (1 108 kg.hm-) 小 于 本 研究 的 结果 。 
与 以 往 在 中 国 北方 的 研究 结果 相 比 ， 云 南 草 地 AGB 
平均 密度 普 逗 大 于 北方 ， 如 马 文 红 等 握 ] 利 用 AGB 与 
遥感 数据 之 间 的 关系 ， 计 算出 中 国 北方 草地 AGB 密 
度 为 790 kg:hm ^; 贾 峰 等 "9 将 19 个 北方 草地 地 面 
定位 监测 点 的 实测 资料 数据 进行 了 分 析 、 归 纳 和 总 
结 ， 得 出 中 国 北方 典型 草原 草地 的 AGB 密度 为 
1 112 kghm-?。 我 国 草地 植被 总 AGB 为 29 232 75 t, 
本 研究 结果 显示 云南 生物 量 范 围 为 1 026.86 万 ~ 
1408.54 万 t, 按照 这 一 估算 范围 ， 计 算出 云南 草地 
生物 量 占 全 国 比 重 介 于 359-489, 在 计算 全 国 草 
地 生物 量 时 不 可 忽略 。 

本 研究 利用 植被 群落 特征 信息 提升 了 复杂 地 形 
条 件 下 草地 生物 量 估算 精度 ,但 在 方法 上 仍 存在 一 
些 不 确定 性 。 首 先 , 采用 的 是 1 : 100 万 草地 分 类 地 
图 ,难免 会 忽略 面积 较 小 的 草地 ,造成 总 的 草地 面 
积 偏 小 ， 使 估算 的 总 生物 量 偏 小 。 其 次 , 在 进行 空间 
尺度 转换 时 可 能 会 产生 误差 ("1 使 用 的 是 250 m 空 
间 分 辨 率 的 MODIS 影像 ， 影 像 算法 容易 将 分 散 的 小 
块 草地 和 大 块 草地 的 边缘 判别 为 非 草地 植被 类 型 的 
地 块 , 由 此 计算 的 NDVI 可 能 与 实地 调研 的 生物 量 
数据 不 匹配 ,在 进行 模拟 时 会 高 估 生 物 量 ; 1:100 
万 草地 类 型 图 与 250. m 分 辨 率 的 遥感 影像 在 进行 空 
间 匹 配 时 ， 地 图 边缘 会 产生 像 元 带 来 的 误差 ， 这 种 
误差 是 难以 避免 且 不 容易 评估 的 。 此 外 ， 只 对 整个 
云南 省 的 草地 数据 进行 统一 建 模 ， 并 未 分 植被 类 型 
建 模 ,没有 充分 考虑 不 同 草地 类 型 相对 独立 的 特点 
和 规律 , 因此 建立 的 模型 具有 一 定 的 不 确定 性 。 如 


果 采 用 更 高 分 辨 率 的 遥感 影像 (比如 Landsat 数据 ) 可 
能 会 减 小 尺度 转换 带 来 误差 , 但 可 能 带 来 新 的 问题 ， 
例如 : 整个 云南 省 共 需 17 幅 Landsat 影像 ， 在 进行 
拼接 时 ,难免 出 现 边界 突变 的 状况 , 而 且 经 常 遇 到 
所 需 时 期 的 影像 有 云 遮盖 的 情况 ， 这 将 导致 在 计算 
植被 指数 时 产生 偏差 。 因 此 ,综合 考虑 后 采用 
MODIS 数据 。 利 用 激光 雷达 数据 (LiDAR) 获 取 的 植 
被 真实 高 度 能 让 估算 结果 更 加 精确 PM" 但 这 类 数据 
尚未 大 范围 应 用 于 草地 植被 ， 从 数据 获取 到 数据 处 
理 都 不 成 熟 扩 。 本 文 只 能 采用 简化 处 理 手段 ， 用 不 
同类 型 草地 植被 的 平均 高 度 代 替 真 实 高 度 ， 做 初步 
的 探索 性 研究 ， 精 确 估 算 大 面积 的 草地 AGB， 需 要 


进一步 探索 如 何 从 激光 雷达 数据 或 直 感 立体 影像 中 
提取 植被 的 垂直 特征 信息 。 
4 结论 


本 研究 综合 遥感 数据 与 地 面 实测 生物 量 数据 ， 
采用 优化 的 估算 模型 分 析 了 云南 省 草地 AGB 及 其 
空间 分 布 格局 , 得 出 以 下 结论 : 

1) 用 植被 特征 信息 (高 度 和 盖 度 ) 确 实 优 化 了 草 
地 AGB 遥感 估算 模型 ， 优 化 后 的 模型 比 之 前 的 拟 合 
优 度 提升 0.289， 使 估算 精度 提高 24.9%， 此 方法 更 
适宜 地 势 复杂 、 和 气候 多 变 、 草 地 植被 类 型 分 布 零散 
的 区 域 。 

2) 云 南 省 2012 一 2014 年 草地 年 均 AGB 介 于 
1 026.86 万 ~1 408.5473 t, 平均 为 1221.11 73 t, 占 全 
国 的 4.1%， 表 明 在 计算 全 国生 物 量 时 不 可 忽视 以 云 
南 为 代表 的 南方 草地 。 云南 省 草地 AGB 密度 总 体 上 
呈现 西部 高 东部 低 、 南 部 高 北部 低 的 空间 分 布 格局 。 

3) 云 南 省 各 市 / 州 之 间 草 地 AGB 密度 介 于 1 130.12~ 
2 116.03 kghm ^; 不 同 草地 类 型 之 间 AGB 密度 相差 较 
A, 由 山地 草 旬 到 暖 性 草丛 再 到 热 性 草丛 有 增 大 的 趋 
势 ， 表现 为 : 山地 草 旬 类 (1 071.73 kg-hm ”< 低地 草 旬 
类 (1 552.45 kg-hm-)< 热 性 灌 草丛 类 (1 579.80 kghm ?)« 
上 暖 性 灌 草 从 类 (1588.12 kg hm ?)« BE f£ S& M 2$ 
(1 771.02 kg'hm-)< 热 性 草丛 类 (2 004.37 kg:-hm )。 
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